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We provide a detailed description of a new symmetry structure of the monomial (Slater) expansion 
coefficients of bosonic (fermionic) fractional quantum Hall states first obtained in Ref. [j], which we 
now extend to spin-singlet states. We show that the Haldane-Rezayi spin-singlet state can be 
obtained without exact diagonalization through a differential equation method that we conjecture 
to be generic to other FQH model states. The symmetry rules in Ref. [l| as well as the ones we obtain 
for the spin singlet states allow us to build approximations of FQH states that exhibit increasing 
overlap with the exact state (as a function of system size) . We show that these overlaps reach unity 
in the thermodynamic limit even though our approximation omits more than half of the Hilbert 
space. We show that the product rule is valid for any FQH state which can be written as an 
expectation value of parafermionic operators. 

PACS numbers: 73.43.1, 11.25.Hf 
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I. INTRODUCTION 

Our understanding of the physics of the fractional 
quantum Hall effect (FQHE) has benefited greatly from 
the existence of model wave functions. Laughlin's trial 
wave function for the v = 1/3 filled FQH state provided 
the first paradigm to understand the emergent behavior 
of interacting electrons in a strong magnetic field^. The 
current understanding of trial wave functions predomi- 
nantly uses the conformal field theory (CFT) connection 
first proposed in Ref. H. For every existing CFT, one can 
build a FQH trial wave function by taking the expecta- 
tion value of branch-cut free primary fields in the CFT. 
The Read-Rezayi (RR) states are a product of this line 
of reasoning^. Both spin-polarized as well as spin-singlet 
states can be obtained this way, most prominent exam- 
ples of which are the Haldane-Rezayi (HR), the Non- 
Abelian spin singlet (NASS), and the Halpcrin states^—. 
FQH trial wave functions are essential to understanding 
the physically important concepts of fractional Abelian 
(in the Laughlin and composite fermion states^—) and 
non-Abelian statistics (in the Moore-Read (MR)^ and 
RR states^). 

The central drawback of the CFT-motivated trial wave 
function approach is the lack of both an explicit decom- 
position of a trial state in a second quantized many- 
body basis and of a first quantized closed form expression 
for the state. As a consequence, Monte Carlo methods, 
while useful for Laughlin states 2 -, cannot be applied for 
most non-Abelian states. Any quantitative analysis of 
these trial states has hence so far relied on exact diag- 
onalization (ED) method s 11 ' 12 . In these methods, one 
starts with a trial Hamiltonian and generates the (low- 
est) Landau level (LLL) Hilbert space. The computa- 
tional effort of diagonalization depends algebraically on 
the Hilbert space dimension, which grows factorially with 
system size. This sets the size limit that is reachable from 



ED. It is hence essential to use all available symmetries 
contained in the trial state and in the associated trial 
Hamiltonian to find the smallest subblock structure of 
the Hamiltonian matrix in terms of the non-interacting 
basis. One symmetry is the reflection of angular mo- 
mentum L z — > — L z which, for a sphere geometry, is 
equivalent to the indistinguishability of the north and 
south pole. Other symmetries such as total L 2 multiplet 
structure exist in some cases. However, they are rather 
obvious in general and do not gain us deep insight in 
the structure of the FQH states. For the Laughlin 1/3 



state, previous attempt; 
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to calculate the weights of 



the free many-body wave functions in the full interacting 
state failed. These worksi^ obtain only 0(1/N\) of the 
0(N\) coefficients, and hence represent a set of measure 
zero of the Laughlin state expansion. In a recent paper, 
two of us found that a large series of FQH trial states 
obey a new type of symmetry for their free basis expan- 
sion coefficients^. The symmetry relates a subset of the 
coefficients of the expansion in free many-body states of 
a given FQH state to products of state coefficients from 
smaller system size. This was developed for bosonic and 
spin-polarized fermionic states. In particular, it was ob- 
served in Ref. 0] that the overlap of the exact FQH state 
with the state approximated by the "product rule" sym- 
metry increases with system size asymptotically towards 
unity. 

In this paper, we give a detailed account of a general 
differential equation method used in Ref. U to access the 
monomial (Slater) decomposition of bosonic (fermionic) 
FQH states. We provide a detailed description of the 
very condensed derivation in Ref. [l] of the expansion coef- 
ficients for bosonic and polarized fermionic states. From 
there, we explain how the trial state can be numerically 
generated at a level intended for the novice reader. Next, 
we present an extended proof of the product rule sym- 
metry for FQH trial states (and for all Jack polynomials) 
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previously summarized in Ref. [l|. We then extend the 
product rule symmetry, which allows to generate even 
more expansion coefficients than previously allowed. We 
also generalize the whole approach to spinful trial states, 
and illustrate it in detail for the HR state. We first 
derive an annihilation operator for the HR state from 
which we develop a recurrence relation for the expansion 
coefficients. Wc investigate the product rule symmetry 
analogue for spinful states and extract the entanglement 
spectrum of the HR state. For a spinful trial state, we 
find that the particle number N, angular momentum L, 
and spin multiplet S are the quantum numbers of the 
reduced density matrix subblocks. 

The article is organized as follows. In Section m wc 
discuss the recurrence relation of Jack polynomials that 
leads to the monomial expansion coefficients for bosonic 
FQH trial states. We elaborate on numerical subtleties 
for certain negative Jack parameters a. In these cases, 
denominator divergences appear in the recurrence for- 
mula: they are accompanied by an (at least) similarly 
vanishing numerator. In Section llHl the Slater expansion 
coefficients of spin-polarized fermionic FQH states are 
derived from a fermionic version of the Laplace-Beltrami 
operator. This is the expanded version of previous calcu- 
lations presented in Ref. OJ. The approach is used to de- 
velop a recurrence formula for fermionic FQH trial states. 
In Section HVt we provide a largely expanded proof of the 
product rule symmetry. For non-Abelian bosonic states, 
we extend the product rule to treat general cases of cut- 
ting through a multiply occupied root partition orbital. 
In Section [Vj we generalize the entire approach to the 
spinful HR state. We derive the recurrence formula, show 
the product rule property, and compute the entangle- 
ment spectrum of this spinful trial state. In Section IVTl 
we take a general viewpoint on the product rule sym- 
metry from conformal field theory. We show that the 
product rule manifests itself as a generic property of all 
FQH states which can be written as an expectation value 
of parafermionic operators hence including a large set of 
both spin polarized and spin unpolarized FQH states. 
Finally, we conjecture in Section IVHI that the product 
rule symmetry is a structural property of the majority 
of FQH trial states including fermionic or bosonic states 
and spin-polarized or spin-unpolarized states, and poten- 
tially serves as an important ingredient to density matrix 
rcnormalization group approaches for FQH systems. 



II. BOSONIC STATES 

FQH states are analytic functions of the positions 
of electrons in a magnetic field. The single-particle 
orbitals in the Landau Level are given by 4> m (z) = 
(27rm!2 m ) -1 / 2 z m exp(— \z\ 2 /4) with angular momentum 
L z = mh, although from now we will neglect the triv- 
ial Gaussian multiplication factors. A non-interacting 
N-particle basis state can be indexed by a partition A 
- an ordered list of the L z angular momentum of the 



10 10 10 1 [6, 4, 2, 0] 
1001110 [6, 3, 2, 1] 

FIG. 1. Pictorial example of the squeezing operation in oc- 
cupation language (left) and partition language (right). The 
squeezing operation takes the first row into the second. 



occupied orbitals. The corresponding occupation num- 
ber configuration is n(A) = {n m (\),m = 0, 1, 2, . . . } 15 i 16 , 
where m labels the individual single-particle orbitals and 
n m {\) is the multiplicity of orbital m in A. We consider 
FQH states decomposed in this many-body basis, either 
of bosons (permanents) or fermions (slaters) with expan- 
sion coefficients c\. The central task of this paper is to 
develop methods to compute these expansion coefficients. 

Wc now define a two-body operation on the many-body 
basis that is important for the purpose of the paper: 
for a pair of particles in the orbitals mi and m-2, with 
mi < itt-2 — 1) the elementary squeezing operation con- 
sists of the two particles shifted to different momentum 
orbitals as n mi 2 -» n mi 2 - I, n TOl 2 ±i -> n mi 2 ±i + I . 
This means that both particles in the mi,m2 orbitals 
are shifted " inwards" the partition (as shown in Fig. [IJ . 
The squeezing defines a partial ordering relation between 
two partitions A > /i when /i is generated by squeezing 
operations acting on A 17 ' 18 . This ordering yields a tree 
hierarchy a complete example of which is shown in Ap- 
pendix [Bj By contrast, when A and /i do not relate by 
squeezing, no ordering relation is set between these par- 
titions. 

The trial FQH states we consider are all squeezed poly- 
nomials. They possess a unique partition, called the root 
partition, dominating all other partitions. This means 
that all partitions with possible (but not guaranteed) 
non-zero weight are generated by subsequent squeezing 
operations acting on the root partition. In many cases, 
this already allows us to omit a significant (more than 
half) part of the Hilbert space (see Table I|. 

In this Section, we focus on the bosonic FQH 
states. The non-interacting basis is given by monomi- 
als M\(zi, ...,z N ) = Per(z l Aj )/n m ^m(A)!, where Xj is 
the momentum index of the jth particle in the partition 
A and Per is the permanent. It was showiJ^ that the 
iV-particle bosonic RR k series of states (which includes 
the Laughlin and MR state) arc a special class of sym- 
metric polynomials. Specifically, this class is called the 
r = 2 single Jack polynomials J"(zi, . . . , Zjv) of parame- 
ter a = and root partition A = [k0 r ~ 1 k . . . k0 r ^ 1 k}. 
The Jack wave functions can be related to WA^_i confor- 
mal field theories and can be classified in terms of sym- 
metric polynomial categories^—. Moreover, the quasi- 
particle excitations of the trial state systems can also 
be written as coherent state superpositions of Jack a 26 ' 27 . 
This provides a complementary view to that of other ap- 
proaches for FQH quasi-particlc excitation states 2 ^—. 
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Algorithmic steps to generate the Jack state 

(i) Generate squeezed monomial basis M K \ Mi — M\ is the root partition 

Order basis states by the integers Yli 2 Ai+JV_I 

(ii) Compute all p K (a) (Eq. [4} 

(iii) Loop over all k: Loop over all pairs of elements li < lj 6 k; 

For each pair unsqueeze to upper dominant partitions p > k and read off c\fj, 

(iv) Compute contribution to c\ K by Eq. [3] 

if Pre (a) = px{o) compute the limit prescription lim € _K}[a — > a — e] 



TABLE II. Sketched howto to use Eq. to generate bosonic FQH states in monomial basis and Jacks of arbitrary parameter 
a. 



nbr particles 


full dim. 


squeezed dim. 


4 


18 


16 


5 


73 


59 


6 


338 


247 


7 


1656 


1111 


8 


8512 


5302 


9 


45207 


26376 


10 


246448 


135670 


11 


1371535 


716542 


12 


7764392 


3868142 


13 


44585180 


21265884 


14 


259140928 


118741369 


15 


1521967986 


671906876 


16 


9020077206 


3846342253 


17 


53885028921 


22243294360 



TABLE I. Size of the monomial basis for the bosonic Laugh- 
lin state v = 1/2 up to N = 17 particles. The second col- 
umn is the complete size. The third column is the number of 
partitions allowed by the squeezing operation from the root 
partition 1010101 . . .0101. 



We expand the Jacks into the monomial basis : 

JS = ^c XK (a)M K , (2) 



K<A 



where K runs over all monomial partitions squeezed from 
the root partition A. There is a known recurrence relation 
for the expansion coefficients cak(q) 50 ' 51 : 



2/ a 



£ {{k+t)-{lj-t))c^{a), 



Px {a) - p K (a) k< ^ x 

(3) 

where k = [l±, . . . , U, . . . , lj, . . . ] and fi = [Zi, . . . , Zj + 
t, . . . ,lj — t, . . . ] denote partitions. We arrange the mo- 
mentum orbitals denoted above in decreasing order from 
left to right, i.e. l\ > I2 > U > lj . . . in k, and a possible 
rearrangement occurs in fi depending on t. All partitions 
(i are understood to be reordered in this way. The sum 
in ([3]) extends over all partitions pi strictly dominating 
k but being dominated (squeezed from) or equal to A 
that can be generated by unsqueezing (i.e. the inverse 
operation to squeezing). The p's are defined as: 



Px{a) = J2 X * (Ai-l--(i-l) 



(4) 



Jacks are eigenstates of the Laplace Beltrami (LB) op- 
erator—: 



1 \ " Zi 

l<3 



+ Za 



' % dzi J dzj 



(1) 

Until recently, the spectrum of the LB operator had been 
studied in detail for only positive a in the context of 
the Calogero-Sutherland models. Recent progress has 
shown that the LB operator has well-defined polynomial 
solutions for certain negative In particular, as stated 
above, some of them are found to correspond to bosonic 
FQH trial states for the ground state, quasi-electron and 
quasi-hole excitations^. 



We now explain an easily implementable computer al- 
gorithm that allows one to obtain a large number of 
bosonic FQH states with high precision. From Q we can 
uniquely index every partition by ^ j; 2 Ai+JV ~\ For any 
numerical implementation, we order the basis according 
to this index, so that the look-up of a partition in the ba- 
sis list becomes logarithmic in effort. By recurrence, we 
can always compute the coefficient of a partition from 
those coefficients of the partitions dominating it. The 
number of dominating partitions is small. Averaged over 
all partitions, it approximately scales with the number 
of fluxes (~ number of orbitals) times the square of the 
number of particles i.e. iV 2 7V0. Thus, to compute all ex- 
pansion coefficients, the procedure gives linear effort in 
the monomial basis dimension. The algorithm to gener- 
ate the Jack state is sketched in Table HU 
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In contrast to positive a for which this formula was 
originally derived, there are minor caveats for certain 
negative a. Situations occur in which the denominator 
in ([3]) vanishes for certain partitions. In these cases, one 
can find that two different partitions /Ltj and /X2 obey 
p Ml (a) = (a), a situation that can be proved to never 
arise for positive a. An elementary example would be 
the 6-particle partitions fJ,i(— 4/3) = [5,5,4,1,1,0] and 
M2(— 4/3) = [3,3,3,3,2,2]. However, this denominator 
divergence is always regularized by a vanishing numera- 
tor. Under a limiting prescription lim^ol ^ a ~ e L the 
quotient either gives or a rational value. Numerically, 
we let a slightly deviate from its exact value, vary it, and 
find a plateau value, which is then identified as the re- 
sulting expansion coefficient for the partition. This type 
of p degeneracy does happen neither for the r = 2 RR 
series nor for the Gaffnian stat o 40 ' 52 ' 53 . It only occurs for 
Jacks that were shown not to be uniquely defined by a 
single clustering condition 4 ^. 

III. POLARIZED FERMIONIC STATES 

Similar to the bosonic case in Sec. UH we start with 
single particle orbitals of the Landau level defined be- 
fore, i.e. <j> m {z) = (27rm!2 m )- 1 / 2 z m exp(~|z| 2 /4) with 
angular momentum L z = rah and m the labeling index 
for all single particle orbitals. However, for the many- 
body state, we now assume that the particles described 
by the first quantized wave functions obey fermionic 
statistics. As a consequence, the non-interacting free 
fermion basis is given by Slater determinant states: sl\ = 
A z (z^z^ 2 ...z^ N ) = Bet(z^). sl x is the unnormalized 
orthogonal Slater determinant, where A denotes the an- 
tisymmetrization over all z coordinates. Different nor- 
malizations can be applied to put the polynomial wave 
function on different manifolds such as the plane or the 
sphere. As in the bosonic case, we describe the free many 
body states by partitions (or occupation numbers). We 
again assume the partition A = [Ai , . . . , Xn] to be ordered 
by decreasing order in angular momentum A; of the ith 
particle. As before, the squeezing operation shifts two 
particles inwards (towards each other) in the partition. 



For fermions, multiple occupancy is forbidden due to the 
Pauli principle. 

In first quantized notation, bosonic and fermionic trial 
states can be transformed into each other by multiplica- 
tion with a Vandermonde determinant. In terms of single 
particle coordinates, this polynomial is the Jastrow fac- 
tor, which is the antisymmetric homogeneous polynomial 
of degree 1. Starting from a Jack polynomial J", the 
transformation reads J" — > S 1 " := J" Oi<j ~~ z j)- The 
polynomials are the exact fermionic analogue of the 
bosonic (Jack) trial state J". For example, the v = 1/2 
bosonic Laughlin state (Jack of (k,r) = (1,2)) becomes 
the v = 1/3 fermionic Laughlin state. The filling al- 
ways changes from bosonic filling v = p/q to fermionic 
filling v = p/{p + q). However, in second quantized no- 
tation, multiplication by the Vandermonde determinant 
does not transform a single monomial to a single Slater. 
To obtain a one-to-one correspondence between a bosonic 
basis and fermionic Slaters, one would first have to trans- 
form the monomials to Schur functions^ 4 -. However, this 
involves knowledge of all the Kostka numbers, a long- 
standing unsolved mathematical problem with no known 
efficient algorithm^. There are two ways to remedy this 
problem. First, we can use the knowledge that the trans- 
formation from monomials to Schur functions is exactly 
given by the J a=l Jack polynomial coefficients, which we 
can compute from Eq. [3l 

However, we try to tackle the fermionic trial states 
from a different angle. We define the decomposition of 
the S" polynomials into Slaters: 

N 

S%( Zl , ...,z N ) = J A « B Y[(zi - Zj ) = J2 (5) 

i<j fi<\ 

To avoid confusion, Ab denotes the bosonic root partition 
and A the associated fermionic root partition. All parti- 
tions y, are squeezed from the fermionic partition A that 
is related to the bosonic partition by A^ = Af + (N — i). 
We now use that the Jack part of is an eigenstate of 
the LB operator, i.e. HlbJ\ b = E\ B (a)J" B . We then 
relate the derivatives acting on J" b to derivatives on S" 
(details are given in Appendix [C|l: 
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E XB (a)S% = l[(z k -z l ) 



k<l 



y> / d \ | 1 Zj + z 3 ( d _ d 
^\ % dzi) a ^ .:, V. * 3 <):., 



E 



Zi 9 1 ^ Zi + Zj f d d 



i,m z<j J 



' dzi dzj 



+ 



E z i( z i H~ z m) 



i,m 



i,m,n 



\%i %m) \ %i z ri) 



1 Zi + Zj I 

a^.Zi — z ■ ' ^ 



2 *J 



c 2 



E 



(0) 



Simplifying several polynomial sums that yield constants as shown in Appendix ^ we can define a fermionic Laplace 
Beltrami operator that diagonalizes 5", i.e. 'Hf jB (a)S"(zi, . . . , z/v) = E\(a)S%(zi, ■ ■ ■ , %n), with 



E x (a) = E A * ( Al_2 (^ 1 ) i ) + (a"" 1 ) ((^^E^ - ^^ -1 )) 



(7) 



1 ( 1 



(8) 



r 



We now diagonalize the above operator in the basis of we realize that, due to its two-body nature, the action of 

Slater determinants. The action of the kinetic part yields Hj on any Slater determinant decomposes into the sum 

HksI^ — ) S V" wnere the [ii denotes the poly- of two-particle interaction terms. It is thus sufficient to 

nomial power of the zth particle in the partition. The ac- look at the action on the two-particle Slater determinant 

tion of the interaction part Hi is non-diagonal in Slater s\ fl= ^ pi ^ = z^ 1 z% 2 — z^ 1 z^ 2 . Assume [i\ > /i 2 , and 

determinant basis and demands further calculation. First define p = — \i<i'- 

I 



ffjSl(p 1)(!l2 ) Z\+Zl,p , m n Z\ + zf 



Z l Z l' 



Zl - 22 



■(*? + *!) -2 



— fp(^ +1 +4 +1 +^2 + ^l 

— Z? V 



zi - z 2 



(21 - z 2 y 

p/2 

)-*£(*? 



1 jp-2 s+1 , p-s+2 a -l , p-s+2 s -l 



1 p/2 / \ 



Zl - z 2 

H J. 

p/2 (s-l)/2 

2^r +2 -4- s+2 ) E (^-*- 1 4- 1 +^- 



p/2 



s/2 



,s-t„t-l 



i )(9) 



s=l t=l s=l t=l 

The two terms are already grouped to yield two-particle Slater determinants. Collecting all prefactors, this gives: 



(Ml-A'2)/2 

fl/sl^,^) = (m -M2 -2)81(^,^3) +2 (/ii-/i 2 -2s)sl ((til _ SjM2+s ). 



(10) 



r 



Eg. 1101 has a particular form: it only scatters "inwards" decreasing relative momentum of the particles, and thus 
the two-particle basis of Slater determinants, i.e. towards to a squeezed partition. Let us now look at the total 
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Occupation Number Angular Momentum N=4 Laughlin State 
9876543210 Q = — 2 



1001001001 



0100110010 



=[9,6,3,0] 



1,5,4,1] 



unsqueezed partitions 6: 



Pa=[9,G,3,0] - 21C 



tion of the partition coefficient for the v = 1/3 Laughlin 
state is shown in Fig. [5]). The complete decomposition 
of the N = 4 particle v = 1/3 Laughlin state is given 
in Appendix [B] As in the bosonic case, this is an ad- 
vance in the numerical computation of the coefficients: 
the computational effort required to generate the state 
scales linearly with basis size. This approach has already 
been applied to increase the maximally reachable system 
in finite size studies^. 



1000110001 V- 2161208 


[(S 


- l)fr[9,5,4,0] 


1000 101010 H 


-(8 


- 4)&[9 i 5,3, 1 ] 


01 010 10001 H 


"(5 


- 1) & [8,6,4,0] 


0101001010 H 


"(5 


- 4)&[ 8 ,6,3,1] 


oiiooooho H 


"(5 


- 4)&[8,7,2,1]] 


FIG. 2. The recurrence relation 


M 


| for the N 



cle v = 1/3 Laughlin state. The partitions are written in 
decreasing order of orbital angular momentum, which ranges 
from 9 to in the case considered. The coefficient of the par- 
tition /j, — 0100110010 is computed with the knowledge of the 
coefficient for the partitions dominating fi. 

action of H^ B on S? expanded in Slaters. The above 
scattering Hamiltonian and the linear independence of 
Slater determinants provide a recurrence relation for the 
coefficients in ([5]) . We collect all diagonal terms and 
invert the sum over s in Eq. [TU] to a sum over all domi- 
nating partitions: 

where pf(a) = X^Ai(Aj + 2i(l — 1/a)). Similar to the 
bosonic recurrence formula in the sum in (jTTj) ex- 
tends over all partitions 9 = . . . , /ij + s, . . . , fij — 
s, . . . , /ijv] that dominate the partition fj, = \p,i, . . . , hn] 
and are squeezed from the root partition A. A new factor 
(— 1)A SW appears: a sign according to the even/oddness 
of the number of transpositions (swaps) of particles from 
a given dominating partition 9 back to fj.. This term 
appears since the reordering of the partition in Slater de- 
terminant language may cause a minus sign due to the 
fermionic anticommutation relations. iVsw starts from 
zero for partition [i and advances by one unit every time 
the momentum of the unsqueezed electron passes through 
the value of the momentum for another electron. As a 
further difference from the bosonic Jack recurrence rela- 
tion, the terms summed in Eq. QT] do not explicitly de- 
pend on the partition 9. This is because the rescaling 
of the s in ([TTJ|) exactly cancels the term's dependence 
on s. For a = — (k + 1), (fTTj) gives the coefficients of 
the fermionic Read-Rezayi states (an example computa- 



IV. PRODUCT RULES 

The coefficients of the monomials (Slaters) in the 
bosonic (fermionic) FQH states exhibit a hidden sym- 
metry found in Ref. [l| and named product rule. The 
product rule is valid in any quantum mechanical nor- 
malization, be it on the plane, sphere, cylinder, disk or 
any other genus-0 geometry. It is however easiest to see 
and explain in the basis of Eq. [3] and Eq. [TTJ for which 
we have already developed the formalism of the previous 
sections. Once a state is obtained in this basis, it is only 
a matter of specific change in the normalization of the 
free many-body wave functions to go between different 
genus-0 geometries. The product rule found ini and ex- 
plained in detail in this section is valid not only for FQH 
Jacks but for all Jacks at any a, of any partition A. Fur- 
thermore, as shown in Section |VT1 below, the product rule 
is a property of quantum Hall trial states even beyond 
Jack polynomials^—. 

We consider a Jack state generated from ^ (or (fTT|) ) 
that can serve as a suitable example to demonstrate 
the product rule. We discuss the fermionic MR state 
of N = 10 particles. This state can be written as 
a linear superposition of Slater determinants squeezed 
from the root partition n(A) = 110011001100110011. 
We pick a configuration squeezed from A that has the 
special property that two parts of the partition can be 
identified as squeezed from individual root partitions 
for smaller systems sizes (see Fig. [5]). Let us consider 
101101000101111010. We observe that the first 7 or- 
bitals from the left, i.e. 1011010, can be squeezed from 
the N = 4 partition 1100110. The remainder right part, 
i.e. 00101111010, can be squeezed from the N = 6 par- 
tition 01100110011. We find that the product of the two 
coefficients obtained from the N = 6 part and the discon- 
nected N = 4 part (—70 and —2 respectively) gives the 
coefficient (-70) • (-2) = 140 of the N = 10 partition. 
The product rule (symmetry) allows the computation of 
a certain set of coefficients of an A-particle state from the 
knowledge of the state for N — 1 particles. This hints at 
similarities with Feynman disconnected diagram summa- 
tion in interacting systems, where the total contribution 
is given by the product of the disconnected components 
(Fig. [3]). As we show below, the product rule approxi- 
mation succeeds in keeping the essential part of the cor- 
relation of the FQH state. We first prove the product 
rule for fermionic Jacks by induction principle; a similar 
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proof can be obtained for the bosonic Jacks. 

Basic induction case. Assume we start with any 
fermionic polynomial S*{zi, . ■ ■ , Zn) of a configuration 
pb < A that can be divided in two disconnected sets: p, = 
(pa,Pb), with Na particles in the first subpart A and 
Nb = N — Na particles in the second subpart B. This 
means that pa is squeezable from an A^-particle root 
partition A^ = [Ai, . . . , \n a ], an d Pb is squeezed from 
the iV^-particle root partition A# = [Ajv a +i, • • • , Ajv]- 
The basic induction case is given by any partition of the 
monomials in S x (z%, . . . , zjy) for which the product rule 
holds. Trivially, for this purpose we can choose the root 
partition A. By definition, it has coefficient 1, and we 
can think of it as being separated into any product of 
two subpart root partitions. Again, all these have coef- 
ficient 1 by definition, so that the product rule holds for 
the root partition itself. 

Induction hypothesis. We now assume the product rule 
is valid for all partitions sharing the separable form 
= (6a, Ob), where pa < Oa < Xa and pb < Ob < A_b. 
As shown in Eq. [TU1 the coefficients of S? are given as 
a recursion from partitions that dominate p. By con- 

I 



struction, any partition dominating p and entering (|10[) 
is also separable according to = (Oa, Ob). As the un- 
squeezed operation is a two-body operation, the sum 
over all dominating partitions 6 can be decomposed into 

£e ;M «KA = J2 t , A <e A <x A +J2„ B <e B <x B: In particular, 
the summation over the individual partition entries pi 
only mixes pi , pj of the left hand side A and right hand 
side B separately, while the remainder right (left) part 
remains unchanged. Partition-wise, the first sum reads 
(9a, Hb), while the second reads (ha, Ob)- Finally, we as- 
sume that all partitions dominating p satisfy the product 
rule: 

bx(e A ,^ B )(a) = b X {0 A ,\ B )(a)b x(XA ^ B) (a), (12) 

where (9a, Xb) denotes the partition formed by Qa and 
the remainder root state partition of part B, Xb- This 
holds vice versa for (A^, Ob). 

Induction proof. Let us consider the coefficient b\^. 
By induction hypothesis, we can rewrite Eq. [TU] as (we 
skip the argument a in the notation for the coefficients 
6): 



= ( E " rf)bx(0 A ,, B} (-l) Nsw + E (»< - ^f) b ^M(-l) Ns ' 



Px - P. 



» \e A ; /i A <e A <\ A 



9b ; Mb <@b < A^ 



2(1-1) 
f IF I °a(a a , Mb ) 



Px - P, 



E (rf-»f)bx(e A ,x B )(-l) Nsw + b 



X(fi A ,X B ) 



E 0*? - rf)h(x A ,e B )(-l) Ns 



9 A ; fj, A <8 A <X A 



9 B ; h b <0 b <X e 



bxp, - -p 



2(i - 1) 



2(i - 1) 



2(i - 1) 



Px (a) - Pl A x b («) + Px (a) - P\ AtlB (a) \ p\ (a) - p^ A Xb (a) p\ (a) - p\ A ^ B (a) 
x E E ~ 4)bf ~ )bx(x A M B) bx { e B ,x A) (-l) Nsw ' A+Nsw - B 

fi A <9 A <X A fj. B <6 B <X B 

bx(^ A ,x B ) ■ bx(x A ,p. B ) q-e.d. 



(13) 



Writing out the coefficients &x(Aa,mb) an< ^ bx(fj. A .x B ) according to Eq. [TU1 we then use p^(ce) — pj t (ce) = p\(c 
Pl A x B (") + Pa (a) - Px A u B («) to S et 



(14) 



A similar line of reasoning applies to the bosonic case. 
We have thus proved the product rule symmetry for this 
type of piece-separable configurations. This is valid for 
all bosonic and fermionic Jack polynomials, and hence 
for all Read-Rezayi states. However, the product rule 
even applies to a much larger range of polynomials (Sec- 
tion IVI[) . Following similar steps as above, the product 
rule can also be explicitly derived for spin-unpolarized 
states such as the Haldane-Rezayi state discussed below 
m Section El 

For bosonic states, there are certain classes of parti- 
tions where it is not immediately clear how to apply 
the product rule symmetry. Let us look again at the 



N = 10 particle MR state in bosonic notation; this state 
is squeezed from the root partition A = 202020202. 

Partitions type I. An easy application of the prod- 
uct exists for configurations such as p\ = 040000600. 
We identify the first 4 orbitals to be squeezed from 
Xa = 2020, while the remainder orbitals are squeezed 
from As = 020202; we find the product rule to hold (see 
also Fig. [3]). 

Partitions type II. Let us analyze the configuration 
Pi = 023000401. We can split the configuration in a 5 
particle separation to the left and a 5 particle separation 
to the right of the cut in Fig. 01 Both parts are discon- 
nected in terms of squeezing operations on the particles. 
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1 1001 lpOl 1001 10011 = 


1100110 


X 


01100110011 


101101000101111010 = 


1011010 


X 


00101111010 


140 


(-2) 


X 


(-70) 


202020202 


202 


X 


020202 


121001410 


121 


X 


001410 


140 


(-2) 


X 


(-70) 






X 





FIG. 3. Product rule applied for a partition of the N — 10 
MR state in fermionic and bosonic notation. The given par- 
tition can be related to the expansion coefficients of product 
of smaller systems, in this case the N = 4 and the N = 6 MR 
state partitions. 



202020202 



023000401 



20201 



02300 



10202 



00401 



J 



-J 



FIG. 4. Generalized product rule for certain partitions of Zk- 
parafermion non-Abelian states in bosonic notation (shown: 
a partition of the N = 10 Moore Read state). In terms of di- 
agrammatic representation of interaction vertices, the orbital 
along the partition cut is assigned a label k for the occupancy 
of the associated root partition orbital. The diagrammatic 
decomposition proceeds by splitting the k particle into j to 
the left and k — j to the right. This way, the product rule 
holds for this class of configurations. 



However, what are the root partitions from which we 
generate the subparts? We have to split one doubly oc- 
cupied orbital of the associated root partition 202020202 
as in Fig. 01 We double copy this orbital and distribute 
the particles in both subparts. We consider 20201 and 
10202 as the root partitions for subpart A and B: to- 
gether they make up 202020202 but the orbital where 
the bold particle 1 is placed is taken to belong to both 
parts A and B. At the same time, we double copy the 
5th orbital of p 2 , i.e. 02300 and 00401 (see Fig. HJ) . Fol- 
lowing this recipe, we find that the product rule holds for 
these type II configurations. We can trace back separable 
configurations of type p 2 to product rule compositions of 
smaller system size (Fig. 0} . 

Let us investigate the accuracy of the approximated 
state for N particles built from the product rule using 
the exact states up to N — 1 particles. To begin with, 



an important quantity is the number of monomials (or 
Slater determinants) whose configurations conform to the 
product rule. Tables IIIII and IIVI show the ratio between 
this number and the total size of the squeezed Hilbert 
space for the Laughlin state (y = 1/2 and v = 1/3) 
and the MR state, respectively. As a rule of thumb, the 
product rule allows to construct more than a third and 
less than a half of the total Hilbert space. This ratio 
decreases with increasing system size, but remains a finite 

> 1/3 fraction of the Hilbert space. 

The important question is how much of the exact state 
is kept in this part of the Hilbert space generated by the 
product rule. We compute the overlap between the ex- 
act state for N particles and the state constructed only 
from the product rule. The overlap is taken using the 
scalar product of the sphere geometry. We find that, 
involving only type I partitions for the Laughlin state 
(Table fll[|) or type I and type II for the MR state (Ta- 
ble EYl, the state approximated by the product rule has 

> 99.9% overlap with the exact state. This tells us that 
the monomials that are generated by the product rule 
contain almost all of the exact state by overlap despite 
comprising only ~ 1/3 of its Hilbert space. In all cases 
we consider, the overlap has the peculiar feature that it 
increases with system size (by contrast, any comparison 
between a model state and the ground state of some re- 
alistic interaction would exhibit the opposite behavior). 
As such, this provides indication that the product rule 
symmetry of quantum Hall trial states becomes exact in 
the thermodynamic limit. Fermionic states show a very 
similar behavior (see Table IIIII for the Laughlin state at 
v = 1/3). The overlaps are also very high but not as 
good as their bosonic counterpart. In our opinion, the 
product rule should be an essential ingredient of future 
DMRG studies. 

For now, we have only considered partitions subject to 
the product rule construction which can be decomposed 
into a pair of ground state partitions of smaller systems 
size. In fact, one can further improve the procedure as we 
find for the specific case of Laughlin states. The prod- 
uct rule can be applied not only when considering dis- 
connected squeezing sequences from the root partition, 
but also from a partition such as 10010001 1| 10001 with a 
cut between the two consecutive particles in the 8th and 
9th orbital. In this case, to reconstruct the Slater de- 
terminant weight, one needs to glue together two Laugh- 
lin states with one quasi-hole excitation each (which are 
also Jack polynomials with the same a parameter as the 
ground state). In this example, this corresponds to con- 
sidering the Jacks of roots 10010001 and 10001. The 
only missing information we additionally need to know 
is the weight of a Slater determinant in a Laughlin state 
which is obtained from the root configuration by a single 
squeezing of two neighboring particles ...1001... into two 
consecutive occupied orbitals ...0110.... It can be shown 
that this weight is always equal to 3 in the basis described 
in Eq. (|10[) . The improvement of the overlap including 
this additional rule is shown in Table [TTTI 
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N 


8 


10 


12 


14 


16 


18 


20 


22 


squeezed dim. 


119 


1070 


10751 


116287 


1326581 


15756587 


193181910 


2429921124 


prod, rules (type I) 
overlap (type I) 


27.73% 
0.8858 


25.23% 
0.9070 


23.58% 
0.9188 


22.48% 
0.9262 


21.63% 
0.9311 


20.95% 
0.9344 


20.40% 
0.9367 


19.95% 
0.9383 


prod, rules (type I + type II) 
overlap (type I + type II) 


72.27% 
0.9875 


70.09% 
0.9895 


68.48% 
0.9919 


66.98% 
0.9930 


65.79% 
0.9936 


64.82% 
0.9941 


64.01% 
0.9944 


63.34% 
0.9946 



TABLE IV. Shown: percentage of Hilbert space that is constructed from the product rule and overlap between the complete 
MR state and the state built from product rule on the sphere geometry. N is the number of particles and the overlap is defined 
as the absolute value of the scalar product. The first row of data is the total dimension of the squeezed Hilbert space for 
different system sizes. The second and third row are the overlap results obtained when only type I partitions are taken into 
account, while fourth and fifth row show the results obtained by involving both type I and type II partitions. 



N 


8 9 10 11 12 13 14 15 16 


v = 1/2 prod, rules 
v = 1/2 overlap 


43.76% 42.01% 40.76% 39.93% 39.52% 39.32% 39.24% 39.19% 39.16% 
0.9933 0.9947 0.9947 0.9956 0.9963 0.9968 0.9972 0.9977 0.9979 


v = 1/3 prod, rules 
v = 1/3 overlap 


47.68% 46.41% 45.33% 44.45% 43.73% 43.11% 42.56% 42.08% 41.65% 
0.9502 0.9534 0.9557 0.9573 0.9585 0.9593 0.9599 0.9603 0.9605 


v = 1/3 prod, rules + qh 
v = 1/3 overlap + qh 


73.92% 72.52% 71.35% 70.40% 69.61% 68.90% 68.27% 67.70% 67.18% 
.9938 0.9944 0.9949 0.9952 0.9955 0.9956 0.9957 0.9958 0.9958 



TABLE III. Shown: Percentage of Hilbert space that can be constructed from the product rule and overlap between the exact 
full state and the state built from the product rule for sphere geometry. N is the number of particles and the overlap is defined 
as the absolute value of the scalar product. The first two rows of data are the results for the v = 1/2 bosonic Laughlin state 
for different numbers of particles. The third and fourth row of data are the results for the v — 1/3 fermionic Laughlin state 
relying only on the knowledge of the Slaters determined by ground state partitions, while the fifth and sixth row additionally 
take into account information stemming from the single quasi-hole excitations. 
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V. HALDANE-REZAYI STATE 
A. Basic properties 

We now turn to the generalization of the bosonic and 
fermionic states involving the spin degree of freedom of 
the constituent particles. In the following, we discuss the 
Haldane-Rezayi (HR) spin singlet stated. The HR state 
was originally proposed as a trial state for the incom- 
pressible plateau state at v = 5/2^. As opposed to the 
spin-polarized Moore-Read (Pfafuan) state at identical 
filling^S, it is a spin singlet. The degree of spin polar- 
ization is still an experimental issue that is not yet settled 
in the v = 5/2 state. The general belief, supported by 
numerical evidence from exact diagonalization studies, is 
that the Moore-Read state is the promising candidate to 
explain the v = 5/2 stated—. The HR state attained 
considerable attention since it can be interpreted as de- 
scribing the transition point between the strong and weak 
pairing phases of a spin-singlet d-wave BCS supercon- 
ductor. This physical interpretation is supported by the 
realization that the bulk CFT is a non- unitary c = — 2 



I 

theory, which is expected to have gapless bulk excita- 
tions^. 

The HR state is given by: 

$hr = A z . 



(zi - wi) 2 (z 2 - w 2 ) 2 ' 

N N N 

x JJ(* - z 3 ) 2 J[( Wi - w 3 ) 2 Y[{zi - wj) 2 , 

i<j i<j i,j 

(15) 

where N is the number of spin f and spin 4- particles 
with positions denoted by z and w, A z , w is the antisym- 
metrizer over the z particles and w particles separately. 
Due to the antisymmetric non-Z7(l) prefactor in the first 
line of (fl~5|) and the evenness of the Jastrow factors, the 
state in (|15[) is fermionic although a bosonic variant can 
also be written. Let us look at the clustering conditions 
at the level of first quantization. The wave function dies 
as the 2nd power of the difference between two equal spin 
coordinates. However, the prefactor of the Jastrow fac- 
tors removes the 2nd order inter-spin zeroes induced by 
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the last Jastrow factor, allowing for configurations where 
one t and one 4- spin sit on the same site. 

In order to span the Hilbert space of spinful states, 
we again start from the single particle Landau or- 
bitals. However, in this case, the single particle or- 
bitals possess an additional spin quantum number a 
taking on the values j" or \,. The free many-body ba- 
sis is given by spinful Slater determinants sU At ,A B i = 

A(zt' Af , ■ ■ .,z X N N ' M ,w^> B \ . . .,wfr Bi ), A 4 , At and 
label the momenta of the ith f spin and ith J, spin parti- 
cles, respectively. 

The squeezing operation is defined as before, and ap- 
plies equally to squeezing between f or ! particles. The 
root partition of the HR state is given by XQOQX . . . OX, 
where X =ti, i.e. X denotes an orbital occupied by both 
spins. This is consistent with the first quantized cluster- 
ing condition and with the filling v = 1 /2 in the highest 
partially populated Landau level. We consider the de- 
composition of the HR state as Slater determinants: 

$hr = cx A t,\ B ±slx A f,\ B i, ( 16 ) 

\At,>>B-l<xooaxooox... 

where the sum extends over all spinful partitions ob- 
tained by squeezing operations on the root partition 
X000X . . . OX. Let us consider the 4-particle HR root 
configuration X000X. It can be expressed in partition 
language as [4 t, 4 1, f, 1} or [4 f, f] x [4 |, 4-]. The 
latter is a factorization in f-spin and J,-spin partitions. It 
differs from the former by a minus sign. Thus, to avoid 
ambiguities of global minus signs due to different ordcr- 



ings of fcrmionic operators, we order all partition entries 
(i.e. t and 4 spin) first by momenta in decreasing or- 
der. For a given momentum, we write the 'f-spin entry 
before the 4-spin entry. Whenever a factorized partition 
notation of f spin and 1 spin momenta appears in the 
following text, it is only for reasons of presentation. The 
ordering of fcrmions should always be interpreted as ex- 
plained above. Moreover, Slater determinants differing 
by an overall spin rotation have equal coefficients. This 
is so since the HR state is a spin singlet and thus spin 
rotationally invariant. 

B. Differential action for the HR state 

From conformal field theory considerations, $o( z , w ) = 
A.U ^-^p ( Z2 -w 2 yi... ) sat i sn cs the following differen- 
tial equation^: 

Id 2 ^ / 1 _L d_\\ 

2 dz ? J^i V & ~ w ^ 2 Zl ~ w i dw i ) J ° 

(17) 

with the same equation for z o w. Following section 
IIII1 we use this equation to obtain an operator for which 
$hr is an eigenstate. We rewrite the derivatives acting 
on $o as derivatives acting on $hr. This amounts to 
taking into account the additional derivative acting on 
Jastrow factors for both spin species. Using the inter- 
mediate steps explained in Appendix |D] we derive the 
following differential equation: 



if / d 



q \ 2 fd d 



' dw. 



z% + Wk I d d 



Wk- 



2 Zi - w k \ dzi dw k , 



4 \ 

VF3 



Z t + Z.j 



' dzi 



' 8Zn 



Wi +W-; I d 
I w, 



£r-^-^r ) ]+4iV(2^-37V + l) 



wj \ % dwi 3 dwj 



HR 



0. (18) 



The above equation contains only 2-body interactions. 
The interaction terms are symmetric with respect to t - 
spin and 4-spin variables. Both the inter and intraspin 
interaction are of Laplace-Beltrami type and are familiar 
from our previous calculations of the polarized fermionic 
states. 



1. Equal spin action 

Let us first compute the action of the terms in (fT8|) 
consisting of equal spin interactions. Once solved for one 
species, e.g. the t spin variables z, this also applies for 
the 4 spin terms. This part of the HR operator is given 



by: 

Zi — Zj \ 1 dzi 3 dzj ) ' ^ ^ 

which wc previously encountered as part of the fcrmionic 
LB operator. However, the term ~ l/(z{ — Zj ) 2 of the LB 
operator is missing. This is an important difference. It 
implies that the operator in Eq. [19] does not map a single 
general Slater determinant of arbitrary degree into an- 
other Slater determinant. The equivalent of (fTO]) cannot 
be written in the current case by using only the operator 
in Eq. [19] Acting with the operator in Eq. [19] on a sin- 
gle Slater determinant usually leads to a fraction. How- 
ever, the special sum of Slater determinants that com- 
prise the HR state does map back into a sum of Slater 
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determinants. This is so since the HR state has two Jas- 
trow terms in equal spins; they cancel any fractions that 
might appear upon the action of (|T9|) on a Slater (see 
Appendix [D]). A single Jastrow factor would not. Hence 
a mapping of a linear superposition of Slaters constrained 
in this way maps, under action of (|19j) . back to the space 
of Slaters. Details of this are given in Appendix [DJ Wc 
hence force the identity: 



£ 



d_ 

'"dzi 



_d_ 



(20) 

where sl M is defined as the Slater determinant of one spin 
species. The partition of the other spin species is omit- 
ted in typing but should be implicated. Since from the 
expression of the HR state we know that the linear com- 
bination of states that form the HR state have a dou- 
ble zero a M s W 2 i ■ • ■ z m) ~ Y\i..j{ z i — z j) 2 i the forced 
identity (|2U|) is clearly true. The coefficients 6 M are then 
given in terms of the coefficients a M . Wc find 

+2a^ l - f i J )Y / M~V Nsw , (21) 
o<v I 



where = [n 1; . . . , ^ - s, m +1 , . . . , fXj + s, . . . , fi N ], and 
Agw again denotes the number of swaps needed to re- 
order the partition. 



2. Different spin action 



In order to compute the term in (|18[) involving action 
on both spin species, we consider one single spinful Slater 
determinant slx A t-,As4 as defined before. Following the 
calculation detailed in Appendix [Dl the spin- rotated par- 
tition sIa^AaI always appears with the same coefficient. 
In this way, we find that the action of the t^l part of 
the operator also maps back to Slaters: 



Ez l + w k ( d d \ 



(22) 



where 0a, 0b are the two partitions of spin f and J,, re- 
spectively. They are obtained by squeezing only oppo- 
site spins in Xa^,Xbi, s parametrizes the changed par- 
tition component obtained from squeezing spin f with 
spin |: A = [\i,a,---,K,a + s,...,X N . A ] and B = 
[Xi,b, ■ ■ ■ , Aj,b— s, . . . , Xn,b]- Thus, the summation takes 
each possible pairwise combination of spin t and spin 4- 
particles, and squeezes them. Ngwcr denotes the reorder- 
ing swaps of spins from species a upon this inter-spin 
squeezing operation. 



C. Recurrence relation 

We are now ready to compute the full action of the 
operator in Eq. (fT8|) on the Haldane-Rczayi state. The 
remainder terms are of non-interacting kinetic type and 



I 

straightforward to compute. By power counting, wc find: 

N ( 1 8 \ 

£ Wfl=4JV(jV-l)$ M . (23) 

The second order derivative terms give 

£(^V $ HR = £ A i<I>HR. (24) 

In the intra-spin term, the sums over differences of A's 
can be computed as 



N 



N 



£(A* - Xj) = £(JV + 1 - 2i)X h (25) 



while there is no similar closed form expression for the 
inter-spin term. Summing up all terms of (|18[) , we find 
the following equation: 
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E c XAt , XBl sh At , XBi Un(W 2 - AN + 1) - i + A L + 2 * A ^ + 2*A 

A A t,A B 4.<X000X000X... V t 

+ 2 E ( A *> CT - Ai ' ff ) 

E C ^t-A B l ( E ^ - ^Aj)sle^,X B i(-l) N - + E - Aflj)slAxt,9B4(-l) JV8W 

AAt,A B i<XOOOXOOOX... \0 A <A A 8 B <\ B 

+ \ E (V"^^.^ 2 -^^ (26) 

In Eq.[2ni we have grouped all diagonal terms to the left and all interactions to the right. Accordingly, the right hand 
side is made of three parts: parts of sums of spin- 4 partitions squeezed from the spin-4 partition Xa, parts of sums 
of spin- 4- partitions squeezed from the spin- 4 partition A b , and parts of sums of spin- 4 and spin- 4- partitions squeezed 
only from particles of different spin from the partition [Xa tj As 4]- We now define the quantities 

PA A T,A B | = ^E( A A ^ + A l ^ + 2 ^,^ + 2^A B , ^ )-i E (V -*«•,*)• ( 27 ) 

We can then immediately show that the root partition obeys pxooox...x = 2N(3N 2 — AN + 1), and thus corresponds 
to the constant terms on the left side of (|26[) . Now, equating the coefficient of every Slater determinant, we derive a 
recurrence relation for the coefficients of the decomposition: 

1 



PXOOOJfOOOX... — Pti A -X,HBi 



E (/*a,< - Maj + 2t)c» AW • (-l) Wsw 

;A t A <e A <toootooot... e A =[^ A ,i,...,/i A)i -i-t,...,^ A ,3-t,...,/i A ,jv] 

+ E 0*B,i - Vb,j + 2t)c^.e Bi ■ (-1)"™ 

e B ;MB<eB<iooo4.oooi... eB=[/*s,i,— .mb,»+*,— .mb.j-*.— i/*b,n] 

+ ^ E - ^ + 2t)c eAt ,flBl(2 - <Wb<W,/^) ' (-1)^+^ , (28) 

e A ,e s; [p A t,MBl]<[e A ,e B ]<[xoooxooox...] / 

where the last sum over the inter-spin terms extends over all partitions 6a, Ob with the property [9a, Ob] = 
[[jUA,l,/M,2, • • • ,VA,i + t,---,LlA,N],[LlB,l,HB,2,---,HB,j - t, . . . , Hb,n]] with p A ,i > ^B.j , or [6 A ,6 B } = 
[[/U,i,/M,2, ■ ■ ■ ,VA,j -t,...,HA,N},lPB,i,HB,2,---,fJ>B,i+t,...,H B ,N]\ with (i B ,i > I^A,j ■ One explicit example of 
computation is presented in Fig. [5j and a complete 4-particle HR Slater decomposition is shown in Appendix |E] 

I 

In analogy to ^ and ([TU| , the Slater decomposition co- HR recurrence relation that hints at further structure in 

efficients of the HR state can be read off from ([28]) . The the HR coefficients beyond the spinful squeezed Slater 

HR state can be generated with linear effort in Hilbcrt basis. The key observation is that the HR state can be 

space dimension. written as££ 



*HR = ]J(z l - Zjf{Wi - WjfP^Zi, Wj ), (29) 



D. Zero-weight partitions 



l< 3 



where P\ is an angular momentum L = polynomial 

The recurrence relations we encountered for the spin squeezed from the root partition AO AO AO . . .AO A. 

polarized bosonic (fermionic) states in Eq. (Eq. [II)]) Interpreted in this way, the spinful 2A— particle HR 

only produced accidental (about 1/N\ of the total Hilbert partitions are generated from the space of partitions of 

space dimension) zero weights for squeezed monomial P\ squeezed from AOAOA . . . AOA times partitions from 

(Slater) configurations. By contrast, we observe an ex- the A-particlc bosonic v = 1/2 Laughlin factor squeezed 

tensive number of zero weight configurations from the from 1010101 . . . 101 for 4 spins and 4 spins, respectively. 
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Occupation Number Angular Momentum N=4 Haldane-Rezayi State 



4 3 2 1 

xooox 



=[4u,4d,0u,0d] 



,,HR 



^= [4,0] T, [4,0] | 



20 



0X0X0 



=[3u,3d,lu,ld] 



JIR 



^=[3,1]T,[3,1]I 



18 



equal spin unsqueezed partitions 

udOdu duOud C[3,1]T,[3,1]| = 20^18 ( 3 ~ 1 + 2 )( c [4.0]f,[3,f]| + C[3,l]t,[4,0u) 
unequal spin unsqueezed partitions 

+|(3 - 3 + 2) (2 - %,1],[3,1]^3,3)(C[4,1]T,[2,1]J. + C[2,l]t,[4,l]|) 

+|(3 - 1 + 2)2(c[ 4j i] T; [3,o]| + C[3,0]T,[4,1]|) 



uOdXO dOuXO 
udOud duOdu 



OXuOd OXdOu 



+|(1 - 1 + 2)(2 - 5[ 3 ,l],[3,l]^l,l)(c[3,2]T,[3,0]l + C[3,0]T,[3,2]i) 



FIG. 5. Explication of the recurrence relation (|28|l for the N = 4 Haldane-Rezayi state. We have chosen to present the explicit 
computation of the coefficient of partition /i = 0X0X0, which involves all conceptional different terms appearing in (I28|l . 



We investigate the 2N particle HR partition 
X000X000X . . . , where we count the momentum in 
terms of polynomial powers beginning from m = at 
the north pole extending to m = 4(A — 1) at the south 
pole. We define the momentum imbalance between total 
t spin and J. spin particle momentum within a partition 
A: | 5^ A* f — Si^i-il- The total momentum summed 
over f and | spins, i.e. Ylio-^hv is 4A(A — 1). For 
N = 2, the partitions of maximum momentum imbal- 
ance are given by J, X f and j 1 j (both with 
momentum imbalance 4 + 2 — (2 + 0) =4). They are 
obtained by just squeezing t spins to the left and I spins 
to the right, and vice versa. If we arrange all t and 4- 
spins to generate the highest momentum imbalance be- 
tween the spin species (which means to arrange all t and 
4- spins around opposite poles) , we find a partition struc- 
ture I 1 ... | 0X0 | | ■ • • t Thus, we find a max- 
imum imbalance of 2N(N — 1) for partitions squeezed 
from the HR root partition. 

Let us now construct the HR state starting from Pi. 
Multiplying a Pi partition of a given momentum im- 
balance with a v = 1/2 Laughlin partition for f and 
4 spins does not change the momentum imbalance, as 
the total momentum added to both spin species is equal. 
Thus, the available range of possible momentum imbal- 



ance in HR partitions is given by the partitions squeezed 
from X0X0X0 . . .X0X. There, arranging the different 
spin species to opposite poles gives a partition structure 
4 4 ■ ■ • 4 X j" j" . . . '|\ The maximum momentum imbal- 
ance is N(N — 1), and thus only one half of the maximum 
imbalance for partitions squeezed from the HR root par- 
tition. As a consequence, we can remove all squeezed 
partitions in the HR state with momentum imbalance 
> N(N — 1), as they must have zero weight. For N = 2, 
this applies for the partition | X | and its spin- 
rotated counterpart | 1 j, as the momentum im- 
balance is 4, while the maximum allowed momentum im- 
balance is 2 (see Appendix [Ej) . While this rule is easily 
implemented numerically, there are even more weight 
partitions squeezed from the HR root partition, which 
cannot be written in a similar closed form as the mo- 
mentum imbalance constraint. Rather, we have to ex- 
plicitly check whether a partition squeezed from the HR 
root partition can be generated from Pi times Laughlin 
partitions. We have sketched the algorithm in Tab. [V] 
by which partitions squeezed from the HR root partition 
can be efficiently denied or accepted. 

E. Spinful Entanglement spectra 



We now focus our attention on the entanglement spec- 



trum, a quantity recently used to identify the topolog- 
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Algorithmic steps to accept / deny a 2N— particle HR-squeezed partition p 
For a given partition p 

(i) Check whether total momentum imbalance | £V fii t f — ih,i\ is < iV(7V — 1). 

(ii) Loop over all f and J. ^V — particle Laughlin partitions k 

(11.1) Divide partition p by «: Check whether all momenta are still in allowed range and no equal spin momenta are mapped 

onto each other (fermionic state) 

(11.2) Check whether the resulting partition can be squeezed from Pi 



TABLE V. Sketched algorithm to remove zero-weight partitions squeezed from the HR root partition XOOOXOOO^l .... We 
apply the structure of the HR state according to (|29[) . By this procedure, the total Hilbert space is reduced to about 83% of 
the full squeezed Hilbert space. 
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FIG. 6. (Color online) Entanglement spectrum for the N = 12 
Haldane-Rezayi state, half sphere cut, Na ~ 6, = sec- 
tor, entanglement levels plotted versus the angular momen- 
tum L^. Levels relating to different spin multiplets, denoted 
by S A , are given in red (singlets) , crossed green (triplets) , and 
square blue (quintuplets). A large part of the main entangle- 
ment weight resides in the singlet sector. 



ical fingerprint of quantum Hall states^—, supercon- 
ductors^, and topological insulators^—, and to detect 
non-local order in spin systems^— and Bose-Einstein 
condensates*^. We cut the sphere into two parts and 
separate the orbitals into regions A and B. From there, 
we compute the reduced density matrix of region A as the 
partial trace of p over the degrees of freedom in region 
B, i.e. pa = Tr^/O. We define entanglement levels £ as 
the spectrum of an entanglement Hamiltonian Ha that 
relates to the reduced density matrix by pa = e~ Hji . De- 
pending on the basis for which the partial trace is defined, 
the entanglement levels possess certain quantum numbers 
associated with commuting operators [p^Qa] = 0. For 
the HR state, the levels are specified by quantum num- 
bers for the number of particles, angular momentum, and 
spin multiplet in A, i.e. Qa — Na, Lf, S A . As the full 
HR state is a spin singlet, it follows that the reduced 
density matrix not only commutes with Sf, but also 

S A such that pa decomposes into spin multiplets. The 
N = 12 HR spectrum for the half sphere cut with half 
the particles in region A is shown in Fig. [51 In case of an 
even number of sites in region A which gives a decomposi- 



tion into integer spin multiplets, we always choose to look 
at the Sf = spin sector where we can study all total 
spin multiplets in one sector. Analyzing the single high- 
est angular momentum entanglement level which also is 
a spin singlet (see Lf = 36 in Fig. [6]), we find that the 
eigenstate of this entanglement level is exactly the N = 6 
HR state. This demonstrates the product rule property 
of the HR state (the fact that the highest angular mo- 
mentum entanglement level of a trial state corresponds 
to the trial state itself for the number of particles in the 
remainder region A is found for all trial states where the 
product rule applies). The product rule can also be ob- 
tained from the explicit recurrence relation between the 
coefficients for the HR state (Eq. |2"5)) by performing sim- 
ilar steps to the ones used in proving the product rule 
for spin-polarized states. The second highest weight sec- 
tor has a single entanglement level that is degenerate in 
entanglement energy with the highest momentum state. 
Again, for a fermionic slater determinant description of 
quantum Hall trial states, this degeneracy is constrained 
by symmetry of how to arrange the slater with one less 
fermionic angular momentum^, and hence can again be 
understood from the Slater determinant structure of the 
HR state. As shown in Fig. [6j we find that both various 
different angular momentum sectors and spin sectors sig- 
nificantly contribute to the dominant entanglement levels 
of the HR state. A significant part of entanglement of the 
HR state is encoded in the spin degrees of freedom: the 
HR state is a spin singlet that cannot decomposed into a 
product of local singlets. 

VI. CONFORMAL FIELD THEORY 

We investigate the product rule property of FQH states 
from the perspective of conformal field theory (CFT). 
For the FQHE, we are concerned with chiral CFTs in 
2 dimensions: all the fields arc holomorphic, i.e. they 
depend on z = x+iy, and not on z. This is an illustration 
of broken time-reversal symmetry in this context. For an 
in-depth introduction to CFT, we engage the reader to 
Ref. [9l|. In this Section we will present the key notions 
that are needed to derive the product rule of quantum 
Hall wave functions from CFT. 

The conformal dimensions arc quantities that appear 
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in the two-point correlation functions of the fields 3>i 
present in the CFT. Conformal invariance fixes any two 
point function to be of the form 



($i(zi)5> 2 (> 2 )) 



Ci, 



(30) 



and its scaling properties are described by the so-called 
conformal dimension Aj. 

Primary fields are of special importance in a CFT. 
Among the infinite number of fields of a given CFT, pri- 
mary fields can be thought of as highest weight fields. All 
non-primary fields (" descendants" ) are created by acting 
on primary fields with certain types of lowering operators 
- the (negative) Virasoro modes. The two point function 
of primary fields $i and $j reads: 



J %3 



{zi - z 2 ) 



2Ai 



(31) 



It is non vanishing only if $ 3 and $j arc conjugate from 
one another, namely if the identity field 1 is produced 
in their fusion. The identity is a very special field, and 
for our purposes it will be the only primary field with 
a vanishing conformal dimension. As a consequence it 
is its own conjugate, 1 = 1. In some cases, such as the 
logarithmic c = —2 ghost system corresponding to the 
Haldanc-Rezayi state encountered in Section [Vj this is 
not the case, and one has to treat this cases more care- 
fully as discussed below. 

Conformal invariance also fixes the form of the corre- 
lation function between 3 primary fields: 



($l(zi)$ 2 (Z2)$3(z3)> = 

-Dl,2,3 



\A-2A 2 



(32) 



( Zl - z 2 )*- 2 *s(z 2 - z 3 ) A - 2Al (^i - z 3 y 

where A = A x + A 2 + A 3 . When the fusion coefficient 
-Di,2,3 is non-zero, the field $3 is said to be produced 
in the fusion $1 x <f> 2 . The fusion rules encode which 
primary fields are produced in any fusion. For instance 
a fusion rule of the form 



$j x $ 2 = $ a + 3> fc + <J> C 



(33) 



tells us that the fusion of $1 and $2 only produces three 
primaries: $ Q ,$b and <I> C . Implicitly, all their descen- 
dants will also appear. 

In order to compute correlation functions, fusion rules 
are not sufficient. An operator product expansion (OPE) 
is a refinement of the fusion rules. It is a formal (and 
exact) expansion of the product of two fields, and is an 
implicit way to define any correlation functions of a CFT 
by iteratively reducing the number of fields involved. It 
has the following generic form: 



$i(zi)$ 2 (z 2 ) = 



(zi - z 2 ) 



A1+A2 — A fc 



-v^), (34) 



where the sum is formally over all fields of the CFT. 
Thanks to conformal invariance, it is sufficient to know 



only the singular terms - those for which Aj+A 3 -— A& > 
- and there will be a finite number of them for any ra- 
tional CFT. TV-point OPEs are a straightforward gener- 
alization: 

$ 1 {z 1 )---$ N (z N )=^2F k (z 1 ,--- ,z N )Q k (0). (35) 



We are now in a position to derive the key property 
of CFT correlation functions that is responsible for the 
product rule. Consider a generic TV = TV4 + TV^-point 
correlation function of primary fields. For simplicity we 
denote the complex coordinates of the fields by Xi and 



C{xi\vj) = 

■ ■ • $ NA {x NA )$N A +i{yi) 



■*n(vn b )). (36) 



A particle cut between Na and Nb particles is directly re- 
lated to the asymptotic behavior of the correlation func- 
tion as we spatially separate the two sets of variables Xi 
and yj. We consider the limit 7 — > 00 of C[xi\yyj). Using 
the TV4 points OPE ([35]) for $i(a;i) ■ ■ ■ $n a (xn a ) around 
we obtain 

C(xi\7Vj) = 

X>fc(a>i)(**(0)$Ar A +i(7Vi) ■ ■ • $jv(7J/jv b )>- ( 37 ) 

k 

Thus conformal invariance tells us that the correlation 

function ($ fe (0)$jv A +i(72/i) ■ • • *w(72/JV B )) - if it is non- 
vanishing - scales as ^~ An a+ 1 Aiv-A fc ag ^ _^ ^ 

From this we can immediately infer that the behavior is 

dominated by the (primary) field <I> a (and its conjugate 

5>a)j which is the field with lowest conformal dimension 

A a appearing in both fusions: 



$1 x • • • x $ Na -> $ Q 
$at 4+ i x • • • x $ w -> $a- 



(38) 
(39) 



In addition, for a primary field $ a the function F a ap- 
pearing in the TV^ points OPE (|35|) is given by the cor- 
relation function 

F a {xi,- ■■ ,X Na ) = ($1(2:1) • • •$Ar 4 (xj V/1 )4 , a(oo)) 



lim x 2A "($i(a;i) • • • ^ Na (x Na )^ & (x)). 



(40) 



Finally we obtain the asymptotic behavior for C(xi\jyj) 
as 7 — > 00 given by 

1 A NA + l + -+A N+ A aC{x ^ lyj) „ 
($l(x 1 )---$> NA (x NA )$a,(oo)} X 

{$ a (0)3> NA+1 ( yi )---<t> N (y NB )) (41) 

for any correlation function of primary fields C(xi\yj) 
(f3"o) . Using first principles of CFT, we hence have derived 
a factorized form of the correlation function spatially sep- 
arating two sets of variables Xi and yj . We will see in the 
following that this property implies both squeezing and 
the product rule for FQH wave functions. 
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A. Spinless FQH wave functions and parafermionic 

CFTs 



1. Parafermionic CFTs 

Many fully polarized QH wave functions can be writ- 
ten in terms of parafermionic CFTs. This includes 
Laughlin 2 -, Moore-Read^ and the Read-Rezayi states^ 
(as well as all Jack states^), but also any generalized 
parafermionic states, such as N — 1 superconformal or 
S3 wave function s] 24 ' 57 . 

Here we constrain ourselves to a condensed derivation 
required to obtain the product rule from CFT. More de- 
tails about parafermions in the context of FQHE can be 
found in Ref. [24|. These parafermionic CFTs, denoted 
as contain a set of k parafermionic primary fields 

{*S>o = 1, • • • ^fc-i} with the following fusion rules: 



n+m 



3d k 



(42) 



The conformal dimension of the field is A„ = 

r n(k-n) ^ aQ( j CQn j u g a ^ e jg lj/ fc _ n . 

Parafermionic FQHE wave functions take the following 
form: 

P( Z1 , ■■■,z N ) = (*(zi) • • • 9(z N )) H(zi - z 3 Y' k . (43) 

Using (|41[) for N = Na + Nb in the case in which all 
the primary fields &i are taken to be parafermionic fields 

(r) 

$ in a Zj. theory, we have: 
yv B A 1+ A aCN7%) „ 

<*(*i) • • • *(s w J*-«(°°)><*«(0)*(l/i) • • • y(y NB )) 

(44) 

where a = Na mod k corresponds to the sector in which 
the Na particles live after the cut, and A a = 1 a(k-a) ^ 
Equivalently, for the wave function this reads P(xi\yj) = 



1 I ' !h ' II- 



P$l(*i)P$l(-Vvj)[Ilvi 



(45) 



In particular, in the neutral sector a = we have 

PN{x l \ 1 y )^^^-)^\[{ 1 y l )^^P NA { X )P NB {y), 

(46) 

where we introduced N^{N) = r(N — k)/k. 

2. Squeezing and product rule 

From (|45[) one can derive the following properties: 
Squeezing. Consider the power of 7 in (|45[) for Ng = 1 : 



This shows that for any monomial entering the de- 
composition of Pn we have /Ui < Ai = f (A 7 — k). By 
iteration on Nb in (|45|) . one finds that (i) the root parti- 
tion is A = (kO r ~ 1 kO r ~ 1 • • • k) and (ii) any other partition 

fi obeys fj,i H h in < Ai + ■ • • + Aj, i.e. is obtained by 

squeezing from A. 

Product rule. In the limit we consider we send the Nb 
particles to infinity: 

lim j-^^^+^^-^^-^PNix^y,) (48) 

7— > 00 

only the monomials m M such that 

Mi H fJ>N B = Ai H \ Nb (49) 

survives. This kills all monomials obtained by squeez- 
ing through the cut between Na and Nb particles, and 
leaves the others invariant. Now that squeezing has been 
cstablished, the product rule is simply equivalent to the 
monomial decomposition of the factorization property: 

lim l- [Xl+ - +XN ^P N {x l \iy ] ) = 

rN A /k 



N B 



4i(^i(%)(n% 



(50) 



where Pn b stands for the north (south) pole reflection of 
Pn b - 



B. Spin singlet states 

The same argument from above applies to spin- 
unpolarizcd states, such as the Haldanc-Rczayi 5 -, 
Halperin£, and NASS state* 7 -. Their CFT description 
involves several types of electron operators, typically con- 
sisting of a parafcrmion field and a vertex operator of a 
set of chiral boson fields. These boson fields $ c and <3? s 
usually describe charge and spin associated with the par- 
ticles. The electron creation operators take the generic 
form: 

y T (2) = * t (2):e75( v ^^ +v ^* s ) :, (51) 
Vi(w) = ViH : e ^(v^*-v^* s ) (52) 

where *S?f and ^ can be trivial fields (Halpcrin) , Gepner 
parafermions (NASS)^ or ghosts (Haldanc-Rezayi)^ 5 -. 
The values of the rational numbers /? and 7 entering the 
vertex operators depend on this CFT. 

The spin polarized wave function assumes the form 

P(zi,Wj) = (* t (zi)* ; (wi) ■ • • * t (zjv)*4(tUiv)) 



ZjY{Wi 



w 



(53) 



1 -j 



Pn(zu z 2 , ■ ■ ■ ZN-1,JZn) 



7 



r(N-k)/k 



(47) 



Taking the asymptotic factorized behavior as n up spin 
and m down spin electrons are taken to infinity as in (|4ip , 
we first obtain a weak form of squeezing: there exists an 
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integer N(n,m) dominating all the partitions (/z^,/^) in 
the sense that 



Mi + ---Ml + 4 + ---/«m <N{n,m), 
which is due to the asymptotic behavior P ~ 7 



(54) 

N(n,m) 



This integer N(n,m) can be expressed in terms of the 
CFT data a, /3, A t , A 4 and A Q . 

Second, we obtain a product rule if a partition /i = 
(/!-)-, /i4.) is separable, i.e. /i = ha + A's- With /ia 
(p-b) squeezed from Xa (Ab), the corresponding mono- 
mial (slater) coefficient is — mj 1(11 x ni\ 2fl2 . In 
order to obtain a stronger form of squeezing, a more de- 
tailed analysis is involved. One would have to specify 
the OPEs of the operators and ^4., and in particular 
work out the dimension of the field $ a appearing in the 
fusion of n fields ty^ and m fields 4^, as we did for the 
spin-polarized case. The detailed analysis of the NASS 
and Halpcrin states is beyond the scope of this paper. 
We treat the HR state in Appendix [F] 
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Appendix A: Helpful Formulas 

We present several formulas used to simplify certain 
polynomial terms. They follow from elementary algebra 
and are easily derived by exploiting the symmetry with 
respect to the summation indices. In Sec. IIII1 we use 



VII. CONCLUSION 

In this paper we have given an extended derivation 
of a recurrence formula for the Slater decomposition of 
fcrmionic Jack polynomial states. We have given a rig- 
orous and detailed account on the product rule symme- 
try for spin-polarized quantum Hall trial states first pre- 
sented in Ref. [l|. We generalized the whole approach to 
spinful states and specifically derived a recurrence rela- 
tion for the spinful slater determinant decomposition of 
the Haldane Rezayi state for which we have computed its 
spinful geometric entanglement spectrum. The product 
rule symmetry is found to be a deep general property 
of quantum Hall trial states, involving both fcrmionic 
and bosonic as well as spinful and spin-polarized states. 
While for states described by parafermionic conformal 
field theory (which include but transcend the Jack poly- 
nomials) we were able to show that the product rule 
comes out of the parafermionic fusion properties, for 
many other FQH states the product rule holds even 
though the states cannot be described by analogue con- 
formal field theory. We showed that the product rule 
can be used as an increasingly good approximation of 
the FQH state. 

While the product rule symmetry is not exactly valid 
for ground states built from realistic interactions (such as 
the Coulomb potential), we would like to investigate if an 
approximated version of this symmetry also manifests it- 
self in these cases. Such property would greatly help to 
improve any density matrix renormalization group algo- 
rithm designed for fractional quantum Hall systems. 
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ZjZj - W k (Zj - Zj) 
{Zi - W k )(zi -Wl) ' (Zi - W k )(Zj - w k ) 

WjWj - Z k {w l + Wj) 

Ll (wi - z k ){wi - zi) ' (w l - z k )(w 3 - zk) 



E 



E 



2N 2 (N - 1). 



(A3) 



Appendix B: Example for Monomial decomposition 

In Fig. [7] we give an example of how to use the squeez- 
ing and product rule properties to obtain the N = 4- 
particle v = 1/3 Laughlin state. 



Appendix C: Fermionic Laplace-Beltrami operator 

We want to rephrase the derivatives on J" s as a deriva- 
tive action on S". This can be done in a compact form 
since they only differ by a Jastrow factor multiplication. 
The first derivative yields: 



z—S a 

dzi 



z—J a 



k<l 



(CI) 
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Similarly, plugging in the previous result for the first 
derivative, the second derivative action can be rewritten 



Appendix D: Haldane Rezayi state 



1. Derivative action 



l d Zl ' J ^ 



~]_(z k - zi) 



k<l 



^ ' 7.1 - - 



dzi 



E 

m 



E 



(C2) 



Starting from the CFT differential equation for the 
non- unitary term (|17[) . we perform a calculation similar 
to the previous derivation of the fermionic LB operator. 
For generality, we can consider the generalized Haldanc- 
Rezayi state: 



4> 



m.m.n 



- $ no* - tor ii^i ~ %) m - *iT- 



X.l 



(Dl) 

Following a tedious, but straightforward calculation, we 
find 



2?( z< a« i ) 2E 



.^-E^+^^r-^-^E 



m 2 — m — 2 



E 



i¥=3 

n 2 -n-2 



. z% Zj dzi 



>E 



E 

k,i 



n + m 



(zi - w k y 



E 



Zi - w k dzi 



E 



Zi - w k dw k 



£ki.i ~ w k)( z i ~ w i) 



2 ( Zl ~ w fe)(^ - 



6 



(D2) 



One observes that the m = 2, n = 2 case is special: the 
three-body terms have an identical coefficient equal to 3. 
These 3-body terms reduce to a constant once the above 
equation for z is added to the similar equation for w. By 
use of (I A3 p . this leads to the expression (fT8"l) . 



2. Recurrence formula 

a. Equal spin action 



equality, we find 

(m -n 2 + 2)) 

m, «2 

/%ni + l^"2 1 ,ni + l.i2'\ _ 
1^1 2 ~ ^~ 2 z l j ~ 

(,"7711,7712 °mi+l,roj-lA*l ^2 ' ^2 *1 /• 
mi ,m2 

(D4) 



From the differential equation (fl"8"|). we know that un- 
der the action of the equal and different-spin terms, the 
sum over Slaters contained in the HR state must yield 
another sum of Slaters. Thus, to deduce a decomposi- 
tion formula we enforce this condition. Without loss of 
generality, we can consider the operator action on a two- 
particle Slater state, and explicitly demand 




711 ,712 



where we use n\ (n 2 ) as momentum indices of the two- 
particle same-spin Slater determinants. Working out the 
derivative action and expanding the expression inside the 



This equation can be worked out iterativcly, starting 
from the maximum polynomial degree. In terms of the 
b coefficients in the final Slater superposition, the coeffi- 
cients a are given by 

&im,fU2 = 0«i -^2)0^1^2 + 2 E a Mi+».^-i(A t i -M2 + 2i). 

(D5) 

Since the interaction is pairwise, this can be directly gen- 
eralized to larger particle numbers of Slaters, and finally 
yields expression (j2"T|) . 

We now illustrate that the degree of the Jastrow fac- 
tor in the polynomial determines whether the action of 
Eq. ID3I causes fractions or not. As an example, we con- 
sider the sh^o = z\ — z\ two-particle Slater determinant. 
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The action of the operator in Eg. ID3I gives: 



gi + z 2 
zi - z 2 











-2 - 



dzi dz 2 



sL-s.o = 3 



Z\ + Z-2 
Z\ - Z 2 



(4 + 4), 



(D6) 

which cannot be decomposed into polynomials without 
fractions. As stated before, this case does not occur for 
the HR state, since there is an equal-spin Jastrow factor 
of 2nd power. Let us now consider the polynomial which 
is constructed from s^o times a Jastrow factor: 



z\ + z 2 



z\ 

Z\ 



22 
22 







8 



Zla ^2-^— Sl 3i0 (2l - Z 2 ) 



' dz 2 



(3(2? 



21 - Z 2 

+4-4)- 



4)(zi 



z 2 ) + z x z 2 {z\ - z\) 



(D7) 



Since all terms from above cancel the zi — z 2 fraction, 
the polynomial composed out of a Slater times a Jastrow 
factor gives no fractions upon the action of the operator 
in Eq. ID3I This holds for any higher power of Jastrow 
factors multiplied with the Slater determinant. 



one from the J, spin partition. If A^ > A& (n > m as 
above) take the prefactor as given in (|D8[) . If A^ < A^, 
change pf"M- an d add it (the spin-rotated Slater must 
have the same coefficient). Produce all allowed squeez- 
ings and add the Slater terms. If one pair {i, k) and some 
squeezed partition produces an n = m state (i.e. a dou- 
bly occupied single particle state), and the remainder "f 
spin and -J, spin partition arc equivalent as well, this factor 
must be divided by two (or rather, in our case, changed 
from 2 to 1). We then obtain (j2"2"|) . 



Appendix E: Example for spinful Slater 
decomposition 

In Fig. [51 we give an example of the HR state decom- 
posed into squeezed partition from the root partition. We 
show the full decomposition of the N — 4 HR state into 
spinful partitions squeezed from X000X. We observe an 
example of a partition squeezed from X000X that has 
zero coefficient, i.e. uOXOd: it cannot be constructed 
from Pi = XOX times 101 for 'f spins and 101 for I 
spins. 



b. Inter-spin action 

To compute the action on the inter-spin term, we again 
constrain ourselves to a two-particle mixed-spin Slater 
partition z™u>™ + and assume n > m (we remem- 

ber that the m > n Slater has the same coefficient): 

Z\ + Wl ( 8 8 \ , n m , m n\ 

— — z " w " + z Twi) 

z\ — wi \ azi owi J 

(m-n)/2-l 

+2(n -m) (4~ l w? +l + zr +i ^i _< ) 

i=l 

+2(n-m)z[ n+m)/2 w^ +m)/2 . (D8) 

For (n — m) = lmod2, the upper limit of the sum is a 
half integer. In this case, the sum is evaluated as an an- 
alytical extension from an integer to a half integer upper 
limit. To each of these Slaters, we also add the coun- 
terpart n O m. To prevent double counting, whenever 
m = n the counterpart must not be added. In the lan- 
guage of single Slaters, this demands the division of the 
equal-momentum Slater prefactor by 2 in comparison to 
other Slaters. For the general particle case, this yields 
the following recipe: Sum over all combinations of pairs 
{i, k) with one coordinate from the t spin partition and 



Appendix F: Haldane-Rezayi, CFT analysis 

As mentioned before in Section [V] the CFT corre- 
sponding to the HR state is a logarithmic theory that 
contains two fields 1 and 1 with vanishing conformal di- 
mension^. This particularity makes the CFT treatment 
more complicated. Fortunately, the factorization can al- 
ready be obtained from the explicit form of the correla- 
tion function: 

(de{zi)dS{wi) ■ ■ ■ d9{z N )dd{w N )I) = 

The r.h.s. trivially obeys the following asymptotic be- 
havior as Nb up-spins and Nb down-spins are taken to 
infinity: 

(d0(zi)d6(wi) ■ ■ ■ dO^N^dO^ZN^i) ■ ■■ 86 (jw n a +n b )) ~ 
7~ 2Ar - (8d(zi) ■ ■ ■ 86(w Na )) (86(z Na+1 ) ■ ■ • 86(w Na+Nb )). 

(F2) 

This yields the product rule (for a cut through an empty 
orbital) using the same steps as in the spin-polarized case. 
Cuts through an occupied orbital can be analysed in a 
similar way, by looking at the asymptotic behavior as Nb 
up-spins and Nb + 1 down-spins are sent to infinity. 
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FIG. 7. Monomial decomposition of N = 4 particle v — 1/3 Laughlin state. The particle positions are denoted by 1. The 
coefficients of the Slater partitions are computed according to (|10p . with a = —2. The arrows denote a squeezing relation from 
the upper to the lower partition. In total, there are 4 squeezing levels till the maximally squeezed partition is reached. 
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FIG. 8. Spinful Slater decomposition of the N = 4 Haldane-Rezayi state. The notation corresponds to X =\\. on the same 
orbital, and u =~\, d =\,. The spinful partitions are written in spin-mixed fashion and ordered with respect to momentum. 
The coefficients are computed corresponding to Eq. 1281 The two appearing zero weight partitions are a manifestation of the 
zero- weight rule according to Q29p . The arrows denote that a given partition can be squeezed from the partition above. There 
are three squeezing levels until the maximally squeezed partitions are reached. 



